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Revivals of quantum correlations in composite open quantum systems are a useful dynamical 
feature against detrimental effects of the environment. Their occurrence is attributed to flows 
of quantum information back and forth from systems to quantum environments. However, 
revivals also show up in models where the environment is classical, thus unable to store 
quantum correlations, and forbids system-environment back-action. This phenomenon opens 
basic issues about its interpretation involving the role of classical environments, memory 
effects, collective effects and system-environment correlations. Moreover, an experimental 
realization of back-action-free quantum revivals has applicative relevance as it leads to 
recover quantum resources without resorting to more demanding structured environments 
and correction procedures. Here we introduce a simple two-qubit model suitable to address 
these issues. We then report an all-optical experiment which simulates the model and per- 
mits us to recover and control, against decoherence, quantum correlations without back- 
action. We finally give an interpretation of the phenomenon by establishing the roles of the 
involved parties. 
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n a multipartite quantum system one can identify different 
I kinds of quantum correlations (such as entanglement and 

genuine quantum correlations), each with characteristics 
useful as resources for quantum information processing 1 ' 2 . The 
exploitation of these quantum correlation resources is however 
jeopardized by the detrimental effects of the environment 
surrounding the quantum system. Quantum correlations 
between noninteracting qubits in independent environments 
have been shown to decay asymptotically ' 3 or even to disappear 
at a finite time 4 ' 5 under the action of Markovian noise. Non- 
Markovian noise, like that arising from structured environments 
or from strong couplings, has instead been shown to lead, during 
the evolution, to revivals of quantum correlations 4,6,7 , permitting 
their partial recovery and thus an extension of their exploitation. 

Revivals of quantum correlations, initially present in a 
quantum system, have been explained in terms of repeated 
correlation exchanges between the qubits and their non- 
Markovian quantum environments because of the back-action 
from the environments to the qubits 4,6,8-11 . Recently, it has been 
shown that revivals of quantum correlations may also occur when 
the environment is classical 2,12-16 and it does not back react on 
the quantum system. Moreover, in such a case, the environment 
cannot either store or share the quantum correlations initially 
present in the quantum system. It thus appears that for revivals to 
occur it is not essential to have a non-Markovian quantum 
environment with back-action. This fact however leads to basic 
issues about the interpretation of back-action-free quantum 
revivals. In particular, it is essential to establish: (1) the role of the 
classical environment in reviving both entanglement and purely 
quantum correlations, for instance if it may act as a control 
system for what operation is applied to the qubits; (2) the role 
of collective effects of the environment on the qubits; (3) the role 
of the memory effects; (4) the role of possible system-environment 
correlations. 

In order to address these issues and explain the mechanisms 
underlying the revivals without back-action, we consider a simple 
model that excludes the presence of collective effects of the 
environment on the qubits and also takes into account 
decoherence on the qubit evolution, so to be in a realistic 
scenario. We also require that the model present a system 
dynamics where the correlations are analytically measurable by 
known quantifiers, and moreover, it is realizable by a neat 



experimental setup that avoids any side effects that can affect the 
expected system dynamics and complicate the interpretation. 
Here we introduce a model of two noninteracting qubits, initially 
entangled, where only one qubit is subject to a random external 
classical field with inhomogeneous broadening in its amplitude. 
This model satisfies the above requirements. We then report the 
results of an all-optical experiment that simulates this model, with 
a random classical field mimicked by quantum degrees of 
freedom, and allows us to observe and control revivals of two- 
qubit quantum correlations without system-environment back- 
action. Moreover, we show a direct connection between the 
occurrence of the revivals and the intrinsic non-Markovianity 
(not due to a quantum environment) of the system evolution. 
Finally, we give an interpretation of the phenomenon showing the 
role of the classical environment in reviving quantum 
correlations. 

Results 

Theoretical model. We consider a system of two noninteracting 
initially entangled qubits, namely a and b, where the qubit a is 
isolated while the qubit b resonantly interacts with a random 
external classical field, whose phase is either q> + — + k/2 or 
(p _ = — n/2 with equal probability 1/2 and whose amplitude has 
a Gaussian distribution. This model is illustrated in Fig. la. The 
interaction Hamiltonian between the qubit b and a classical field 
E with phase cp, in the rotating frame at the qubit-field frequency, 
is H— ih(Q/2) (<J + e~ t(p — a _e l(p ) where the coupling constant 
(Rabi frequency) Q is proportional to the field amplitude and a + 
are the raising and lowering operators in the basis {|0 > ,| 1 > }. The 
time evolution operator U(t) — e ~ iHt/ ^ has the matrix form 15 

m _ / cos(Qt/2) e-'*sin(Qf/2)\ m 
cpK) V-^sin(Of/2) cos(Qf/2) )' {) 

Due to the randomness of the field phase acting on the qubit b 
and being the qubit a isolated, the global dynamical map 
of the model of Fig. la is K t p% = (1/2) 

(H fl <8> U (p (t))p2^ a <g> ul(t))> where p% is the initial two- 
qubit state. A relevant property of this map is that it moves inside 
the class of Bell-diagonal states (that is, mixtures of the four Bell 
states), for which analytic expressions of the most used quantifiers 
of correlations are known 2,17,18 . Therefore, we shall choose 
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Figure 1 | Model and experimental setup, (a) Illustration of the model of a random external classical field acting only on the qubit b, whereas qubit a is 
isolated. The random dephaser shifts of n, with probability 1/2, the phase (n/2) of the input field, (b), Experimental setup simulating the physical model of 
panel a. The state preparation process is started by separating the ultraviolet pulses into two paths with the Wollaston prism. The relative amplitude 
between the two paths can be controlled by a HWP1. HWP2 and HWP3, with angles set to 22.5°, are used to control the polarization of the pump light by 
rotating |H> and |V> to (| H>+ | V))/y/2 and (| H) — \ V))/y/2, respectively. Photon pairs are emitted into modes a and b. The birefringence in the beta- 
barium-borate crystals is compensated by the quartz plates (CP). The photon in mode b is then coupled by a single-mode fiber and directed to the 
environment part. The Soleil-Babinet compensator (SBC) and quartz plates (QPs) add a relative phase </>, proportional to their length, between |H> and 
\V}. HWP4 and HWP7, with angles set to 0° introduce a relative n phase between horizontal and vertical polarizations. The angles of HWP5 and HWP6 
are set to be 22.5°, which rotates |H> and |V> to \(\ H)+ | V)) and 4- (| H) - | V)), respectively. Quantum state tomography is implemented with the 
polarizations of the final states analysed by a quarter-wave plate (QWP), a HWP and a polarization beam splitter (PBS) in each arm. The photons are 
detected by single photon avalanche detectors (D1 and D2) with 3nm interference filters (IFs) in front of them. 
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special Bell- diagonal states as initial states of the system. At this 
stage the evolution is cyclic and does not exhibit decoherence. In 
order to introduce decoherence in the system evolution, we 
consider the field suffering a noise source due to signal 
inhomogeneous broadening whose effect is to induce a 
Gaussian distribution in the field amplitude and thus in the 
Rabi oscillation frequency Q of the qubit evolution, that is 
f(Q) = {2/oy/n) exp[ — 4(Q — £l 0 ) 2 /<7 2 } where Q 0 an d <J are the 
central Rabi frequency and the Rabi frequency width, respectively. 
The effect of noise is therefore transferred to the intrinsic 
evolution of the quantum system. The final evolved state p°f is 
obtained by tracing out the Rabi frequency degree of freedom 
from A t p%. 

Experimental setup. The experimental setup simulating the 
model introduced above is displayed in Fig. lb. The quantum 
system is given by two polarized photons, each one being a qubit 
with basis states |H> (horizontal polarization) and \ V} (vertical 
polarization). The preparation part of the setup is able to generate 
polarization- entangled photon pairs in a given Bell- diagonal state 
p™ b . The photon in mode a is directly sent to the state tomo- 
graphy part whereas the photon in mode b is directed to the 
environment part and successively to the state tomography part. 

The environment is given both by the two photon paths 
separated by a beam-splitter, the reflected path p + and the 
transmitted path p _ , and by the measurement process that traces 
over the time difference between the two paths themselves, 
creating a statistical mixture of them with equal probabilities 
(1/2). The effect of the two paths is to make the photon basis 
states |H>, |V> undergo, apart from an unimportant global 
phase factor, the unitary transformations (see second section of 
Methods) 

| H) ^ cos(</>/2) | H) ± /sin(0/2) | V), (2) 

| V) ^±fsin(0/2) | H)+cos(</>/2) | V), 

where 0 = cl>t is the phase difference between \H} and |V> 
introduced by the Soleil-Babinet compensator (SBC) and the 
quartz plates (QPs). co is the photon frequency and x = LAn/c is 
the time difference to cross the optical element (SBC or QP), with 
L being the thickness of the optical element, c the vacuum speed 
of light, An the difference between the indices of refraction of the 
horizontal and vertical polarizations. The action of each path in 
equation (2) corresponds to a rotation U p+ (4>) of the basis states 
{|H>,|V>} of b. The two paths p ± of equation (2) act on the 
photon b exactly as the two time evolution operators U (PT (t) of 
equation (1), respectively, with the correspondences |0> <-> |H>, 
|1> <— ► \Vy and (j) = cot <— ► Qt (which implies Q <— ► co). The 
environment part of the setup simulates the random phase of 
the field acting on qubit b (random dephaser) by the action of 
the beam-splitter that separates the two photon paths, corresponding 
to the effect of the field with either phase, plus the measurement 
process that does not distinguish the two paths p + . This mimics 
the random external field acting on qubit b, giving the overall 

output state Ap£ = (1/2) £ p=p± (1„ ® U p )p*(t a ® 17+ ). Up 
to now we have considered a photon with frequency co. In fact, in 
our experiment the photon has an intrinsic Gaussian frequency 
distribution f(co) = (2/ay/n) exp[ - 4(cd - co 0 ) 2 /a 2 ] (ref. 19), 
where co 0 is the center frequency and a the frequency width. 
This frequency degree of freedom introduces decoherence by a 
decoherence parameter k — ff(co)cos[(j)(co)]dco = cos(co 0 t) 
exp( — <7 2 t 2 /16) (ref. 3). This decoherence source is of the same 
kind as that due to the Rabi frequency distribution in the 
model introduced above. The effective evolved state p°f is then 



obtained by tracing out the photon frequency degree of freedom 
from Ap™. 

We prepare initial Bell diagonal states of the form (see first 
section of Methods) 

p£ = *. + \® + )(® + I +*._ I *-><<!>- I, (3) 

where X + +X_ = 1 and | Q> ± ) = (l/>/2)(| HH)± | VV}). X + 
is arbitrarily determined by the relative power between the two 
pump paths created by the Wollaston prism. As a consequence, 
the two-photon output state is the Bell-diagonal state (see second 
section of Methods) 

P3f =(i + ic)/2(a. + |<d + x<d + | +x_ |o-)(o- |) 

+ (1-k)/2(X + \ V + I I I), 1 j 

where | x ¥ ± ) = (l/V^XI HV) ± | VH». When L is small so 
that <jt <C 1 (exp( — o 2 x /16) ~1), the decoherence parameter 
becomes k& cos(0 = co 0 t) and we have a cyclic and thus 
coherent evolution: this is the case when only the SBC, that allows 
us to adjust small lengths, is used. When the QPs are used the 
decoherent evolution comes out. 

Experimental dynamics of correlations. We investigate the 
dynamics of the various kinds of correlations present in the two- 
qubit system. We quantify the total correlations by the quantum 
mutual information 2 i(p°f) an d the entanglement by 
concurrence 1 C(p°jf) = max{0, T}, where T is a function of 
Pab- Q uantum correlations are instead quantified by the quantum 
discord 2 Q(p° a f) = l{p 0 £)-C{p 0 ab)> where C (Pab) indicates the 
classical correlations (see third section of Methods). 

The results for the coherent evolution, realized by the SBC 
device with the QPs removed, are plotted as a function of the 
relative phase cb — co 0 t in Fig. 2a, starting from p™ b with X + — 0.9. 
Total correlations oscillate with period n. Plateaus and sudden 
changes in the dynamics of classical and quantum correlations 
occur several times, showing their robustness against variations of 
the relative phase. Quantum correlations revive after vanishing at 
about k/2 and 3n/2 (in a 2n period). Entanglement exhibits dark 
periods around the same points and then revives. 

In Fig. 2b, we further show the experimental results for the 
decoherent evolution, where both SBC and QPs are used. We plot 
the envelope dynamics of correlations as a function of the QP 
length L for the same initial state p™ b with ^ + =0.9, by 
measuring the maximum amplitudes of revival of the quantifiers. 
A sudden change in behaviour in the dynamics of the maximum 
values of classical and quantum correlations is found when 
L^82^ 0 (^ 0 = 800nm is the center photon wavelength). The 
maximum of entanglement revivals decreases monotonously and 
completely vanishes at L « 258^ 0 - The theoretical curves exhibit- 
ing these decaying revivals are displayed in the inset of Fig. 2b. 
The coherent evolution in Fig. 2a is part of the decoherent 
evolution in Fig. 2b, which is denoted by the red box in Fig. 2b. 
Although the decay behavior of the maximum values of 
correlations is similar to that in ref. 3, the internal noise 
mechanism is completely different. There are no revivals of 
correlations in the Markovian evolution of ref. 3, whereas here the 
correlations exhibit collapses and revivals during the evolution in 
absence of system-environment back-action. 

Our experimental setup controls X + and thus the revival 
amplitude. Figure 3 shows the results for the coherent evolution 
at a relative phase (j) — 7i/4 (also applicable at (j) — 3k 14 in the 
region of first revival). It is seen that all correlations are 
symmetric to the point ^ + =0.5, which represents a classical 
initial state (zero quantum correlations). Interestingly, there is a 
range of values of 1+ (0 < X + <0.14 and 0.86<X + < 1) where 
quantum correlations remain unchanged at their maximum 
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Figure 2 | Dynamics of correlations in the absence of system-environment back-action. The initial input state p[Ji has A + =0.9. The black solid line, 
dashed red line, dot-dashed blue line and dotted dark cyan line represent, respectively, the theoretical predictions of total 1, classical C, quantum 
correlations Q and concurrence C (the value is set to 0 when T<0, which is represented by the solid dark cyan line). The black, red, blue and dark 
cyan dots represent the corresponding experimental results, (a) The dynamics of correlations in the coherent evolution. Only the SBC is used in the 
experimental setup and the relative phase (in units of n) is (t> = co 0 T. Revivals of entanglement and of quantum correlations are clearly noticed, 
(b) The maximum of the correlations is plotted as a function of the QP length L (in units of /l 0 = 800nm, that is the center photon wavelength). The 
coherent evolution, including oscillations and revivals, in a is a part of the decoherent evolution in b, which is denoted by the red box in b. The inset of b 
shows the theoretical dynamics of the various kinds of correlations with a calculation precision set to 0.U o . Error bars are due to the counting statistics. 
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Figure 3 | Control of the amount of correlations. Different kinds of 
correlations as a function of X + (initial state coefficient) in the case of 
coherent evolution, at a fixed value of relative phase </> = a> 0 T set to n/4. 
Due to the periodic behaviour of the correlations as a function of </> 
observed in Fig. 2a, the results are immediately applicable, for example, at 
the point </> = 37r/4 that is in the region of first revival. The black solid line, 
dashed red line, dot-dashed blue line and dotted dark cyan line represent, 
respectively, the theoretical predictions of total X, classical C, quantum 
correlations Q and concurrence C (the value is set to 0 when T<0, which 
is represented by the solid dark cyan line). Black, red, blue and dark cyan 
dots represent the corresponding experimental results. Error bars are due to 
the counting statistics. There is a good agreement among theoretical curves 
and experimental points. 



value, implying that at (j) = 3k 14 they cannot take a value larger 
than that observed in Fig. 2a, for X + — 0.9. This is connected to a 
peculiar behaviour of the sudden change points of Q. In contrast, 
the degree of entanglement revival at (j) — 3n/4 can be increased 
with respect to that of Fig. 2a, for example, by increasing X + . 

Connection to non-Markovianity. The observation of the revi- 
vals without back-action in our setup can be connected to the 
intrinsic non-Markovianity of the system evolution. In our sys- 
tem, photon b experiences the noisy environment and photon a 
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Figure 4 | Revivals and non-Markovianity. The entanglement evolution, 
quantified by the concurrence C, of the initial maximally entangled state 
| <t>+) = (1/v / 2)( | HH) + | W» for the coherent case as a function of 
the relative phase <fi = co 0 T (in units of n). Red dots represent the 
experimental results of C and the black solid line represents the 
corresponding theoretical prediction C = |cos0|. The orange solid line 
represents the evolution of the trace distance rate of change that is 
positive (signature of non-Markovianity) just when the entanglement 
revivals occur. The dashed dark yellow line represents the boundary 0. 



can be treated as an isolated ancilla. Non-Markovianity then 
immediately follows from the non-monotonous behaviour of 
entanglement evolution obtained starting from the input Bell 
state | 0+> = (l/V^XI HH) + | VV)) (ref. 20), displayed in 
Fig. 4 for the coherent case. This aspect is emphasized if we 
consider the non-Markovianity measure based on the trace 
distance D(t) 8 . Any temporary increase of D(t) is a signature of 
non-Markovianity that can be identified by the regions where 
the rate of change a(t) = dD(t)/dt is positive. For our local 
map of random external fields, within the coherent case where 
(j) — co Q T is the dimensionless time variable, we have 21 
a((j)) = — sgn(cos0)sin0 where sgn(x) is the sign function 
(see last section of Methods). We find <j(4>)> 0 just when 
entanglement revivals occur, as displayed in Fig. 4. Interestingly, 
we also notice that the larger entangling rates of the system 
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correspond to the larger signatures of memory effects (larger 
values of o '(</>)). 

The trace-distance measure indicates the presence of non- 
Markovianity when the distinguishability of the pair of states 
increases with time: this can be simply interpreted as a flow of 
information from the environment back to the system which 
enhances the possibility of distinguishing the two states. We 
however remark that this back-flow of information does not 
necessarily imply back-action and correlation exchange. The 
presence of system- environment back-action means that 
the system affects its dynamics via the environment. When the 
environment dynamics is unaffected by the system this back- 
action on the system cannot occur. This is the case considered in 
our model and in the corresponding experimental setup. 



Interpretation. The connection between the revivals and the 
intrinsic non-Markovianity of the system evolution does not by 
itself provide an explanation of the observed phenomenon but 
only gives a necessary condition for it to occur. A possible 
interpretation of the phenomenon of quantum correlation 
recovery in absence of system- environment back-action can be 
provided when one describes the initial global system-environ- 
ment state p se (0) = p s (0) 0 p e of our setup as a quantum-clas- 
sical state 15 , where the quantum part is the two-photon system 
p 5 (0) = p™ b and the classical part is the environment e acting 
locally on the photon b. The environment can only be in a 
separable mixture of its basis states that in our setup can be 

expressed as p e = Ep= P± ( 1 / 2 ) \p) e (P l> where IP± are the 
basis states corresponding to the two random paths (field phases 
in the theoretical model) of equation (2). By this formalism, we 
are able to quantify the correlations present in the global 
quantum-classical system a-b-e. We also show that the 
environment e and the photon b never become correlated 
during the evolution and that this aspect however does not 
jeopardize the role of the classical environment in reviving the 
quantum correlations initially present in the quantum system. 

For the sake of simplicity let us consider the case when 
decoherence is negligible (analogous argumentations hold even if 
decoherence is present). The two-qubit output state is 
obtained by tracing the overall evolved state p se {(j)) over the 

environment e, where p se (4>) = (H fl <8> U be )p se (0)^ a 0 u\ e ) and 
Ube = J2 P = P± u pW ® | p) e (p | is the unitary operation acting 
on qubit b and its environment e. This way, the same dynamics 
observed above for a-b correlations exhibiting the revivals is 
retrieved. During the time evolution, due to H a 0 I7j, e , the states 
of the classical environment are invariant, the qubit b does not 
affect the environment and the back-action by the environment 
on the qubit is absent. Moreover, a classical environment cannot 
store any quantum correlations on its own. One could then think 
that the presence of revivals is due to some correlation established 
between b and e. The subsystem b-e evolves under the local 
unitary operation U he so that the correlations between b-e are 
invariant. If one traces p se (4>) over the isolated qubit a, it is 
straightforward to show that the qubit b and its environment e 
never become quantum correlated. For example, for an initial 
Bell-diagonal state between a-b like that of equation (3), it results 
that the reduced state of b-e during the evolution is the 
uncorrected state (H&/2) 0 (H e /2). Therefore, the qubit- 
environment correlations do not have any role in the occurrence 
of the phenomenon under investigation. 

From the structure of U be it is clear that the environment only 
has the role of a control system for what unitary operation is 
applied to the qubit b. By memory effects, the environment e 
keeps a classical record of what unitary operation has been 



applied to the qubit b and this happens even without back-action. 
It is the lack of this classical information that makes quantum 
correlations disappear at a given time and it is the recovery of this 
information that makes quantum correlations then revive. The 
information the environment holds about a system is therefore 
due to what action the environment performs on the system. To 
better clarify this aspect, let us consider the evolution in Fig. 4. 
When quantum correlations are zero, for example at (j) = n/2, 
there is a complete lack of classical information about what 
unitary operation the environment is applying to the qubit 
(statistical mixing of the two different unitary operations U p + ); 
successively, when the initial quantum correlations are recovered 
at (j) — 7i, this information becomes known because the two 
unitaries both act as a o x 'bit-flip' operation (they are just the 
same operation). 

Discussion 

In this article we have introduced a simple two-qubit model 
particularly suitable to investigate the issue of the occurrence of 
revivals of quantum correlations despite the absence of back- 
action. In this model only one qubit of the pair locally interacts 
with a classical environment, in the form of a random external 
field with inhomogeneous broadening. The model forbids any 
mechanism both of system-environment back-action and of 
quantum correlation storage. Moreover, this model presents an 
overall decoherent evolution where any collective effect of 
external fields on the two-qubit system is excluded. We have 
then reported the results of a novel all-optical experimental setup 
that simulates the introduced model and permits us to observe 
revivals of quantum correlations (entanglement, quantum 
discord), with a decreasing revival amplitude in time. The 
amplitude of revivals can be controlled by adjusting the 
parameter of the initially prepared state. 

We have also experimentally and theoretically shown, as a 
necessary condition for revivals of quantum correlations to occur, 
the link between the revivals themselves and the intrinsic non- 
Markovian feature of the system evolution. We have provided a 
possible interpretation of the phenomenon in terms of quantum- 
classical system, showing the role had by the classical environ- 
ment as a control system for what action it has on the quantum 
system and then in reviving quantum correlations. The 
introduced model and the corresponding experimental setup 
have thus permitted us to identify an interpretation of the 
phenomenon and to exclude side effects. Our findings clearly 
reveal that the revivals of quantum correlations are a dynamical 
feature of composite open systems independent of the nature 
(classical or quantum) of the environment. 

Finally, our results introduce the possibility to recover and 
control, against decoherence, quantum correlation resources even 
in absence of back-action, without resorting to quantum 
structured environments 22-24 or correction procedures 4 ' 25-27 
and moreover open the way to further studies on the issue of 
quantum correlation revivals in classical environments 15 ' 28 . 

Methods 

State preparation process. As depicted in Fig. lb, ultraviolet pulses are frequency 
doubled from a mode-locked Ti:sapphire laser centered at 800 nm with 130 fs pulse 
width and 76 MHz repetition rate, and pass through a half- wave plate (HWP) 1 and 
a Wollaston prism to be separated into two optical pump paths. The HWP1 rotates 
the polarization of the ultraviolet pulses which are separated by the Wollaston 
prism into two paths. The Wollaston prism transmits horizontal polarization |H> 
and reflects vertical polarization | V) . Therefore, the relative amplitude between the 
two pump paths can be controlled by the HWP1. The HWP2 and HWP3 are used 
to control the polarization of the corresponding pump light. The two pump pulses 
are then combined by a polarization-independent beam-splitter and pump two 
identically cut type-I beta-barium-borate crystals, with their optic axes aligned in 
mutually perpendicular planes 29 , to generate entangled photon pairs. After 
compensating the birefringence effect between \H) and | V> in beta-barium-borate 
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crystals with quartz plates (CP), maximally entangled photon pairs are generated in 
the state | 0 + > = (l/v^XI HH) + | W)) ab by the transmitted pump pulses and 
in the state | O" > = (l/v2)(| HH) - \ VV)) ab by the reflected pump pulses. In 
our experiment, the time difference between the two pump paths is larger than 
the coherence length of the pump light so that the path information is traced out at 
the end of measurement 3 . The prepared state is then the Bell-diagonal state pJJ 
of equation (3). 

Action of the environment part of the setup and output state. If the photon b is 
in the horizontal polarization \H) and enters the environment part of the setup 
(see Fig. lb), it undergoes all the transformations due to each optical component in 

HWP5 

the transmitted path (p _ ) and reflected part (p + ) that give, respectively, | H) — > 
(1/V5)(|H>+ I V)) SB - QP (1/V2)(|H>+^ | V» H - 6 cos(0/2) | H)-i 
MM I V> and | H) ™T | H) ^ (1/V5)(| H) + | V» "^(l/vS) 
(| H)+e^ | V» H ™ cos(0/2) \H)-i sin(0/2) | V) ™ P7 cos(</>/2) | H) + 
i sin(</>/2) | V), where we have omitted the unimportant global phase factor e 1 ^ 12 . 
These are the two unitary transformations for |H> given in equation (2), with 
4> = cot. Analogously, one obtains those for the polarization | V> . In particular, we 

HWP5 

find that the transmitted path transformation is | V) -»• (l/\/2)(| H) - \ V)) 
SBC " QP (l/v^)(| H)-e^ | V» H ™ -isin(0/2) |H)+cos(0/2) | V), whereas 

i n i i r I T r \ HWP4 . __. HWP5 , , r-s 

the reflected path transtormation is | V) — > \ V) — > (1/V2) 

(| H) - | V)) SBC -^ QP (l/v^)(| H) - e<* | V> H ™ P6 - i sin(<£/2) | H) + 

HWP7 

cos(</>/2) | V) — >• f sin((/)/2) | H) + cos($/2) | V). These are the two unitary 
transformations for \V} displayed in equation (2). 

We now consider the evolution of the two-photon state | 0 + ) = (l/y/2) 
(| HH) + | VV)) with a single frequency co. The state at the end of the transmitted 
part and of the reflected part is obtained by applying the transformations of 
equation (2) only on the second qubit b and it reads, respectively, 



| 0 + )^ cos(0/2) | 0 + > ± sin(0/2) | v F + >, 



(5) 



where | V F + > = (l/v^XI HV) + | VH». After tracing out the path information 
(the detection in our experiment traces over the time difference between the two 
paths), the input state |0 + ) then becomes the output mixed state. 

|0+)(0 + | "S* cos 2 (0/2) | 0 + >(0 + | +sin 2 (0/2) |^ + >(^ + | (6) 
Similarly, when the input state is |0~ > the output state becomes 

| <D~><<D- | ^ cos 2 (0/2) | <D~><<D- | +sin 2 (0/2) |^->(^F- |, (7) 

where | x ¥~ ) = (l/y/2)(\ HV) - | VH». Therefore, starting from the effectively 
prepared state pJJ of equation (3), the two-photon output state is 

p7 =cos 2 («/>/2) (X+ |0> + }<<i>+ I +x_ |0>-><<i>- I) 

+ sin 2 (</>/2) (X + |4- + )(T + | |T")(T- |), (8) 

The previous output state corresponds to the case of coherent evolution, where 
only the SBC is used (4> = co 0 t). Considering the Gaussian distribution f(co) of the 
frequency co, the decoherence parameter k = ff(co)cos(j)(co)dco arises by tracing out 
the photon frequency degrees of freedom. As cos 2 (0/2) = (1 + cos0)/2 and 
sin 2 (0/2) = (1 — cos0)/2, the introduction of the decoherence parameter k turns 
the previous state p°jf into the output state p° a f given in equation (4). 

Correlation quantifiers. Total correlations are calculated by the quantum mutual 

information T(p ab ) = S(p a ) + S(p b ) — S(p ab ), where p a and p b are the reduced 
density matrices of p ab whereas S(p) = — tr(p log 2 p) is the von Neumann 
entropy 2 . Classical correlations are defined as the maximal information gained 
about p ab by a measurement on one of the qubits, that is 2 



C(Pal) 



(9) 



where il| is the set of projective measurements performed on qubit b, qj = 
tr ab (UjP ab U^) and pi = tr b (Uj p ab U^) /qj is the post-measurement state of a after 
obtaining the outcome j on b. Quantum correlations (quantum discord) are then 
given by Q(p ab ) = 1{p a b) ~ C(Pab)- Two-qubit entanglement is quantified by 
concurrence 1 C{p ab ) = max{0,r}, where T = ^ - jl 2 ~ ^3 ~ y/)U> and Xj are 
the eigenvalues in decreasing order of the matrix p ab (o y 0 Cy)p* ah (o y <S> o y ) with 
o y denoting the second Pauli matrix and p* ab corresponding to the complex 
conjugate of p ab in the canonical basis {\HH) ,\HV} ,| VH> ,| VV) }. The theoretical 
dynamics of the correlation quantifiers is found by using their known analytical 
expressions for Bell-diagonal states 2 ' 17 . The experimental results are acquired from 
the reconstructed output density matrices (by state tomography). Total correlations 
and concurrence are directly calculated for the output states. For obtaining classical 
and quantum correlations the maximization of equation (9) is required, performed 
by scanning the projective operators {Hj} represented by {cosoc|H> + Aina|V>, 
e~ l ®smu\H} — cosoc|V>} on photon b. The scan precision of oc and 9 is set 

to 71/100. 



Quantifier of non-Markovianity. We consider a maximally entangled input state 
| ® + ) ab = (l/\/2)(| HH) + | VV)), with the photon b subject to the environment 
and photon a isolated. The output state, in the coherent evolution case (0 = co 0 t is 
the dimensionless time variable), becomes 

p° a f = l/2(l + cos0) | 0+>(0+ | + 1/2(1 -cos0) | \|/ + )(\|/ + |. The entangle- 
ment of the final state is quantified by the concurrence C = 2 max 
[1/2(1 + cos0), 1/2(1 - cos0)] - 1, that is C(0) = |cos0|. 

We quantify the memory effects present in the evolution of a system using the 
criterion of non-Markovianity based on the distinguishability of two evolving 
quantum states, as measured by the trace distance D(p 1 (t),p 2 (t)) = (l/2)|pi(t) — 

p 2 {t)\, where | A |= TrVA^A, p f (f) = Afp f (0) (i = 1,2) and A f is a dynamical map 
of the open system 8 . The rate of change of the trace distance is 



<7(0 = dD(Pi(f),P 2 W)/df, 



(10) 



Given two initial states of the system, any temporary increase of the trace distance 
(a(t) > 0) is a signature of memory effects (non-Markovianity of the dynamical 
map of the system). In order to obtain the degree of non-Markovianity one then 
takes the maximum of a(t) over the possible initial states Pi )2 (0) and calculates 8 



JV(A) = 



X ft, 2 (0) 



L>o<Mt). 



Let us now consider the dynamical map of random external fields acting on a 
qubit in the coherent case (small lengths of the QPs, SBC compensator) 

A t p(0) = (1/2) E P =± n / 2 ^cp{t)p{0)ul (t), where is defined in equation (1) 
and Clt = coqT = 0. For two arbitrary initial states, p k (0) (k = 1,2) in the basis {|0 > 
= |H>,|1> =|V>} 



P*(0) : 



Pk 

l-OLk 



(11) 



the trace distance evolves as D(p x (0), p 2 (0)) = {(cc\ — a 2 ) 2 cos 2 0+ | jS x — jS 2 | 2 
[cos 4 (0/2) + sin 4 (0/2)] + [(fil - Plf + (p, - p 2 ) 2 } cos 2 (0/2) sin 2 (0/2)} 



1/2 



It 



follows that the maximum of the rate of change is found for the initial states having 
ocj = a 2 = 1/2 and Pi = - P 2 = i/2, that is the pure states | \|/ 1)2 ) = (l/V2)(\ H) 
± i | V)). In this case the trace distance is D(p!(0),p 2 (0)) = |cos0| = C(0) and the 
corresponding rate of change is cr(0) = — sgn(cos0)sin0 21 . 
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